WEAK APPROXIMATION OF A FRACTIONAL SDE 
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Abstract. In this note, a diffusion approximation result is shown for stochastic dif- 
ferential equations driven by a (Liouville) fractional Brownian motion B with Hurst 
parameter H G (1/3, 1/2). More precisely, we resort to the Kac-Stroock type approxi- 
mation using a Poisson process studied in [H [7] , and our method of proof relies on the 
algebraic integration theory introduced by Gubinelli in fTSl . 



After a decade of efforts [2], [6], [131 EDI EH EHl [27] , it can arguably be said that the basis 
of the stochastic integration theory with respect to a rough path in general, and with 
respect to a fractional Brownian motion (fBm) in particular, has been now settled in a 
rather simple and secure way. This allows in particular to define rigorously and solve 
equations on an arbitrary interval [0,T] with T > 0, of the form: 



stands for a ci-dimensional fBm with Hurst parameter H > 1/4. A question which arises 
naturally in this context is then to try to establish some of the basic properties of the 
process y defined by ([T]), and this global program has already been started as far as 
moments estimates [T3], large deviations [121 [23], or properties of the law [5], [2S] are 
concerned (let us mention at this point that the forthcoming book [H] will give a detailed 
account on most of these topics). 

In the current note, we wish to address another natural problem related to the frac- 
tional diffusion process y defined by ([T]). Indeed, in the case where B is an ordinary 
Brownian motion, one of the most popular method in order to simulate y is the following: 
approximate 5 by a sequence of smooth or piecewise linear functions, say (X^)e>o, which 
converges in law to B, e.g. an interpolated and rescaled random walk. Then see if the 
process y'^ solution of equation ([T]) driven by X'^ converges in law, as a process, to y. This 
kind of result, usually known as diffusion approximation, has been thoroughly studied in 
the literature (see e.g. [ISl [3D1 [SI] ) , since it also shows that equations like ([ID may emerge 
as the limit of a noisy equation driven by a fast oscillating function. The diffusion ap- 
proximation program has also been taken up in the fBm case by Marty in [22] , with some 
random wave problems in mind, but only in the cases where H > 1/2 or the dimension 
d of the fBm is 1. Also note that, in a more general context, strong and weak approxi- 
mations to Gaussian rough paths have been studied systematically by Friz and Victoir in 
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1. Introduction 
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[TU] . Among other results, the following is proved in this latter reference: let (X^)£>o be 
a sequence of (i- dimensional centered Gaussian processes with independent components 
and covariance function W . Let X be another d-dimensional centered Gaussian processes 
with independent components and covariance function R. Assume that all those processes 
admit a rough path of order 2, that converges pointwise to i?, and that R^ is suitably 
dominated in variation norm for some p G [1,2). Then the rough path associated to 
also converges weakly, in 2p-variation norm, to the rough path associated to X. 

This result does not close the diffusion approximation problem for solutions of SDEs 
like ([T]). Indeed, for computational and implementation reasons, the most typical processes 
taken as approximations to B are non Gaussian, and more specifically, are usually based 
on random walks [UllSllEH] or Kac-Stroock's type [H [3, [171 [29] approximations. However, 
the issue of diffusion approximations in a non-Gaussian context has hardly been addressed 
in the literature, and we are only aware of the aforementioned reference |22], as well as 
the recent preprint [H] (which deals with Donsker's theorem in the rough path topology) 
for significant results on the topic. The current article proposes then a natural step 
in this direction, and studies diffusion approximations to ([T]) based on Kac-Stroock's 
approximation to white noise. 

Let us be more specific about the kind of result we will obtain. First of all, we consider 
in the sequel the so-called (i-dimensional Liouville fBm 5, with Hurst parameter H G 
(1/3,1/2), as the driving process of equation ([T]). This is convenient for computational 
reasons (especially for the bounds we use on integration kernels), and is harmless in terms 
of generality, since the difference between the usual fBm and Liouville's one is a finite 
variation process (as shown in [3]). More precisely, we assume that B can be written as 
B = {B^, . . . , B'^), where the i?*'s are d independent centered Gaussian processes of the 
form 

Bi= [ (t - r)^-^rfW;, 
Jo 

for a (i-dimensional Wiener process W = {W^, . . . , W^). As an approximating sequence 
of B, we shall choose (X'^)e>o, where X"^'* is defined as follows, for i = 1, . . . , ci: 

X''%t)= [ {t + e-r)^-^e''\r)dr, (2) 
Jo 

where 

e''\r) = (3) 

for A^*, i = 1, . . . , (i, some independent standard Poisson processes. Let us then consider 
the process y"^ solution to equation ([1]) driven by X^, namely: 

dyl = a{yt)dXt + h{yl)dt, = a G M", tG[0,T]. (4) 

Then our main result is as follows: 

Theorem 1.1. Let {y'^)e>o be the family of processes defined by and let 1/3 < j < H , 
where H is the Hurst parameter of B. Then, as e — > 0, converges in law to the process 
y obtained as the solution to where the convergence takes place in the Holder space 
C^([0,T];R"). 



Observe that we have only considered the case H > 1/3 in the last result. This is of 
course for computational and notational sake, but it should also be mentioned that some 
of our kernel estimates, needed for the convergence in law, heavily rely on the assumption 

> 1/3. On the other hand, the case H > 1/2 follows easily from the results contained 
in [7j, and the case H = 1/2 is precisely Stroock's result [29]. This is why our future 
computations focus on the case 1/3 < H < 1/2. 

The general strategy we shall follow in order to get our main result is rather natural in 
the rough path context: it is a well-known fact that the solution ?/ to ([1]) is a continuous 
function of B and of the Levy area of B (which will be called B^), considered as elements 
of some suitable Holder (or p- variation) spaces. Hence, in order to obtain the convergence 
y^^y in law, it will be sufficient to check the convergence of the corresponding approxi- 
mations and X^'^ in their respective Holder spaces (observe however that X^'*^ is not 
needed, in principle, for the definition of y^). Then the two main technical problems we 
will have to solve are the following: 

(1) First of all, we shall use the simplified version of the rough path formalism, called 
algebraic integration, introduced by Gubinelli in [13] , which will be summarized in 
the next section. In the particular context of weak approximations, this allows us 
to deal with approximations of B and directly, without recurring to discretized 
paths as in [6]. However, the algebraic integration formalism relies on some space 
C^, where k stands for a number of variables in [0,T], and 7 for a Holder type 
exponent. Thus, an important step will be to find a suitable tightness criterion in 
these spaces. For this point, we refer to Section HI 

(2) The convergence of finite dimensional distributions ("fdd" in the sequel) for the 
Levy area B^ will be proved in Section [Sj and will be based on some sharp estimates 
concerning the Kac-Stroock kernel ([3]) that are performed in SectionEl Indeed, this 
latter section is mostly devoted to quantify the distance between f{u)9^{u)du 

and f{u)dWu for a smooth enough function /, in the sense of characteristic 
functions. This constitutes a generalization of [7], which is interesting in its own 
right. 

Here is how our paper is structured: in Section [21 we shall recall the main notions of 
the algebraic integration theory. Then Section [3] will be devoted to the weak convergence, 
divided into the tightness result (Section Hj) and the fdd convergence (Section [5]). Finally, 
Section [6] contains the technical lemmas of the paper. 

2. Background on algebraic integration and fractional SDEs 

This section contains a summary of the algebraic integration introduced in [13], which 
was also used in [2H1 in order to solve and analyze fractional SDEs. We recall its main 
features here, since our approximation result will also be obtained in this setting. 

Let X be a Holder continuous M'^-valued function of order 7, with 1/3 < 7 < 1/2, and 
a : M" — s> M"^'^, b : — > M" be two bounded and smooth functions. We shall consider in 
the sequel the n- dimensional equation 

dyt = a{yt)dxt + b{yt)dt, yo = aeW, tG[0,T]. (5) 

In order to define rigorously and solve this equation, we will need some algebraic and 
analytic notions which are introduced in the next subsection. 



2.1. Increments. We first present the basic algebraic structures which will allow us to 
define a pathwise integral with respect to irregular functions. For an arbitrary real number 
T > 0, a vector space V and an integer > 1 we denote by Ck{V) the set of functions 
g : [0, T]^ — > V such that gtl■■■t^, = whenever tj = tj+i for some i < k — 1. Such a function 
will be called a (/c — 1) -increment, and we will set C^,{V) = Uk>iCk(y). An important 
elementary operator is defined by 

fc+i 

6 : C,{V) - Ck+,{V), {Sg),,..,,^, = Y.(-^)'^'9t,-u-t,^.^ (6) 

1=1 

where ti means that this particular argument is omitted. A fundamental property of S, 
which is easily verified, is that S6 = 0, where 6S is considered as an operator from Ck{V) 
to Ck+2{V). We will denote ZCk{V) = Ck{V) n Ker5 and BCk{V) = Ck{V) n lm6. 

Some simple examples of actions of 6 are obtained for g G CiiV) and h G C2(y). Then, 
for any s,u,t E [0, T], we have 

{Sg)st = gt- gs, and {6h)sut = hst - - Kt- (7) 

Furthermore, it is easily checked that ZCk+iiV) = BCkiV) for any k > 1. In particular, 
the following basic property holds: 

Lemma 2.1. Let k > 1 and h G ZCk+i{V). Then there exists a (non unique) f G CkiV) 
such that h = 6f . 

Observe that Lemma [2?T] implies that all elements h G C2(y) with 6h = can be written 
as h = 6f for some (non unique) / G CiiV). Thus we get a heuristic interpretation of 
S\c2(v)'- it measures how much a given 1-increment is far from being an exact increment 
of a function, i.e., a finite difference. 

Note that our further discussion will mainly rely on /c-increments with k < 2. For the 
simplicity of the exposition, we will assume from now that V = M"^. We measure the 
size of these increments by Holder norms, which are defined in the following way: for 
/ e C2{V) let 



sup and C^{V) = {feC2{Vy,\\f\\,<oo} 



s,te[o,T] 



t - s i" 



Obviously, the usual Holder spaces C^iV) are determined in the following way: for a 
continuous function g G Ci{V) simply set 

\\9\\,= \\S9\\,, (8) 

and we will say that g G Ci{V) iff \\g\\fj, is finite. Note that || ■ ||^ is only a semi-norm on 
Ci{V), but we will work in general on spaces of the type 

CU^) = {g ■■ [0,T] ^V;go = a, \\g\\, < oo} , (9) 

for a given a G V^, on which is a norm. For h G C^^V) set in the same way 

ll^ll7,P = sup ^t'u^ . In (10) 



s,u,te[0,T] 

\h 



\u - spit - u\P 



inf 1 ll^i|lp.,M-ft; h = '^hi,0<p^<fi^, 



where the infimum is taken over all sequences {hi G C^lV)} such that h = hi and for 
all choices of the numbers pi e (0,/x). Then || ■ ||^ is easily seen to be a norm on CsiV), 
and we set 

ci;{v) := {heCsivy, \\h\\,<^}. 

Eventually, let Cl'^{V) = Ufj_yiCi^{V), and note that the same kind of norms can be 
considered on the spaces ZC^^V), leading to the definition of the spaces ZCt^iV) and 

With these notations in mind, the crucial point in the current approach to pathwise 
integration of irregular paths is that the operator 6 can be inverted under mild smoothness 
assumptions. This inverse is called A. The proof of the following proposition may be found 
in [13], and in a more elementary form in [H]: 

Proposition 2.2. There exists a unique linear map A : ZCl+{V) Cl+{V) such that 

6 A = Id^f^i+iy^ and A6 = Id^i+iyy 

In other words, for any h G C\'^{y) such that 6h = there exists a unique g = A{h) G 
C2^{V) such that 5g = h. Furthermore, for any p > 1, the map A is continuous from 
ZCf^iy) to C2 (V^) and we have 

\W\,<i^\\h\\,. hezct;{v). (11) 

Moreover, A has a nice interpretation in terms of generalized Young integrals: 

Corollary 2.3. For any 1-increment g G C2{V) such that Sg G Cl^{V) set 6f = {Id — 
AS)g. Then 

n 

where the limit is over any partition Ugt = {to = s, . . . ,tn = t} of [s, t], whose mesh tends 
to zero. Thus, the 1-increment 6 f is the indefinite integral of the 1-increment g . 



2.2. Weakly controlled paths. This subsection is devoted to the definition of general- 
ized integrals with respect to a rough path of order 2, and to the resolution of equation ([5]). 
Notice that, in the sequel of our paper, we will use both the notations fdg or J'stif dg) 
for the integral of a function / with respect to a given increment dg on the interval [s, t]. 
The second notation J'stif dg) ^ill be used to avoid some cumbersome notations in our 
computations. Observe also that the drift term h is generally harmless if one wants to 
solve the equation ([5]). See e.g. Remark 3.14 in [25]. Hence, we will simply deal with an 
equation of the form 

dyt = cr (yt) dxt, t G [0, T], with yo = a (12) 
in the remainder of this section. 

Before going into the technical details, let us make some heuristic considerations about 
the properties that a solution of equation ([5]) should have. Set at = cr{yt), and suppose 



that y is a solution of ffT^ . with y E for a given 1/3 < k < 7. Then the integral form 
of our equation can be written as 

yt = a+ / audxu, t e [0,T]. (13) 



^0 

Our approach to generalized integrals induces us to work with increments of the form 
{^y)st = Ut — ys instead of (fT^ . However, it is easily checked that one can decompose ffT^ 
into ^ ^ 

{6y)st= / audxu = d-s{Sx)st + Pst, with pst = / {au - d-s)dxu, 

J s J s 

if our integral is linear. We thus have obtained a decomposition of y of the form 6y = 
a5x + p. Let us see, still at a heuristic level, which regularity we can expect for a and r. 
If cr is a C^-function, we have that & is bounded and 

\c^t - crs\ < ||Vcr||oo||l/|U|t - si'', 

where denotes the Holder norm of y defined by ([H]). Hence we have that a belongs 
to and is bounded. As far as p is concerned, it should inherit both the regularities of 
Sa and x, provided that the integral J^icTu — as)dxu = Jg{Sa)sudXu is well defined. Thus, 
one should expect that p E Cf**, and even p E C^^^ ■ To summarize, we have found that a 
solution 5y of the equation should be decomposable into 

5y = aSx + p, with a E bounded and p E C'^'^. (14) 

This is precisely the structure we will demand for a possible solution of (fT2l) : 

Definition 2.4. Let z he a path in (M^) with k < 7 and 2k + 7 > 1. We say that z 
is a controlled path based on x, if zq = a, which is a given initial condition in M^, and 
6z E C2 (M'^) can be decomposed into 

6z = C6x + r, i.e. {6z)st = Cs{Sx)st + Pst, s,t E [0,T], (15) 

with C E (R'=^"') and p is a regular part belonging to C|'*(M'^). The space of controlled 
paths will be denoted by Qn^ai^^), one? a path z E QK,a{^'') should be considered in fact 
as a couple {z,Q- The natural semi-norm on Qn^ai^'') is given by 

with Afig-X'^iV)] defined by ^) and ^/[C, C'l{V)] = supo<,<T |C|y ■ 

Having defined our algebraic and analytic framework, we now can give a sketch of the 
strategy used in [13] in order to solve equation ( fT2l) : 

(1) Verify the stability of Qn^aiR-'') under a smooth map : ^ M". 

(2) Define rigorously the integral J Zudxu = J{zdx) for a controlled path z and com- 
puted its decomposition (|T^ . 

(3) Solve equation (fT^ in the space Q^^ai^^) by a fixed point argument. 
Actually, for the second point one has to assume a priori the following hypothesis on the 
driving rough path, which is standard in rough path type considerations: 

Hypothesis 2.5. The W^-valued '^-Holder path x admits a Levy area, that is a process 
x2 = J{dxdx) E Cl^iW^""^) satisfying 

6x^ = 6x^6x, i.e. [{6yi^)sut] = [Sx']su[Sx^]ut, s,u,t E [0,T], i,j E {1, . . . ,d}. 



Then the following result is proved in [13], using the strategy sketched above: 

Theorem 2.6. Let x be a process satisfying Hypothesis \2.5\ and o" : M" ^ M"-^*^ be a 
function, which is bounded together with its derivatives. Then 

(1) Equation admits a unique solution y in Q^^ai^"') for any k < 7 such that 
2k + 7 > 1. 

(2) The mapping (a, a;,x^) ^ y is continuous from x C^iW^) x Cl^^W^""^) to 

Q.,a(M"). 

We shall see in the next subsection that this general theorem can be applied in the fBm 
context. 

2.3. Application to the fBm. Let B = {B^, . . . , B'^) be a d-dimensional Liouville fBm 
of Hurst index H G (|, |), that is B^, . . . , are d independent centered Gaussian pro- 
cesses of the form 

Bi= I {t-r)^-^dWi., 
Jo 

where W = (W^, . . . , W^) is a (i-dimensional Wiener process. The next lemma will be 
useful all along the paper. 

Lemma 2.7. There exists a positive constant c, depending only on H , such that 

E\Bl - Bl\^ = [it - r)^-5 - (s - r)^-^]'dr + f\t - rf'-^dr < c\t - (16) 
for all t > s > 0. 

Proof. Indeed, it suffices to observe that 

f [{t-rf~^-{s-rf~^Ydr = 
Jo 



< 

and that j'it - rf^-^r = ^^f^. 

□ 

Let S be the set of step-functions on [0,T] with values in R*^. Consider the Hilbert 
space Ti. defined as the closure of S with respect to the scalar product induced by 

d 

((l[o,ii], • • • , l[o,td])' (l[o,si]! • • • 5 l[o,Sd]))-H = Si,ti e [0,T], i = 1, . . . ,d, 

i=l 

where R{t,s) := E[BIBI]. Then a natural representation of the inner product in Ti is 
given via the operator J^, defined from S to L^([0,T]), by: 

J^^it) = {T- tf-'^^it) - - ^) - - t)^-l dr, 




and it can be checked that can be extended as an isometry between 7i and the Hilbert 
space L^([0, T]; M"^). Thus the inner product in can be defined as: 

(V', = (^V', ^ ^) L2([0,T];R<*) • 

The mapping (l[o,ti], • • • , l[o,td]) ^ Sf=i c^^n also be extended into an isometry between 
"H. and the first Gaussian chaos Hi{B) associated with B = {B^, . . . , B'^). We denote this 
isometry hj ip ^ B{ip), and B{ip) is called the Wiener-Ito integral of ip. It is shown in [HI 
page 284] that C7(IR'^) C 7i whenever 7 > 1/2 — if, which allows to define B{ip) for such 
kind of functions. 

We are now ready to prove that Theorem 12.61 can be applied to the Liouville fBm, which 
amounts to check Hypothesis 12.51 

Proposition 2.8. Let B be a d- dimensional Liouville fBm, and suppose that its Hurst 
parameter satisfies H G (1/3, 1/2). Then almost all sample paths of B satisfy Hypothesis 
\2.5\ with any Holder exponent 1/3 < 7 < iJ, and a Levy area given by 

dBu® dB,, I.e. Blit,j)= dBi dBl z, j G {1, . . . , 4, 

J s J s J s 

for < s < t < T . Here, the stochastic integrals are defined as Wiener-Ito integrals when 
i j , while, when i = j , they are simply given by 

ft 



^ dB:£dBl = ^{Bl-Bl)' 



Proof. First of all, it is a classical fact that B G CJ(R'^) for any 1/3 < 7 < i/, when 
S is a Liouville fBm with if > 1/3 (indeed, combine the Kolmogorov-Centsov theorem 
with Lemma 12.71) . Furthermore, we have already mentioned that CJ{W^) C Ti. for any 
7 > 1/2 — iJ. In particular, if ii > 7 > 1/3, the condition 7 > 1/2 — iJ is satisfied 
and, conditionaly to B^ , J^dBl^f^dBl is well-defined for i j, as a Wiener-Ito integral 
with respect to B\ of the form B\(f) for a well-chosen ip. Hence, is almost surely a 
well-defined element of C2(M'^^''). 

Now, simple algebraic computations immediately yield that (5B^ = 6B ® SB. Further- 
more, Lemma [6.41 yields 

E[\BU^,J)\^]<c\t-sr. 

Invoking this inequality and thanks to the fact that is a process in the second chaos 
oi B, on which all norms {p > 1) are equivalent, we get that 

E[\Blit,j)f]<c,\t-s\'^". 

This allows to conclude, thanks to an elaboration of Garsia's lemma which can be found 
in [HI Lemma 4] (and will be recalled at that B^ g C^^iR'^'"^) for any 7 < 1/3. 

This ends the proof. 

□ 



With all these results in hand, we have obtained a reasonable definition of diffusion 
processes driven by a fBm, and we can now proceed to their approximation in law. 



3. Approximating sequence 



In this section, we will introduce our smooth approximation of B, namely X^, which 
shall converge in law to B. This will allow to interpret equation (jll) in the usual Lebesgue- 
Stieltjes sense. We will then study the convergence in law of the process solution to 
(jlj) towards the solution y oi 

As mentioned in the introduction, the approximation of B we shall deal with is defined 
as follows, for z = 1, . . . , d: 

X'''{t)= [ {t + e - r)^-h'''{r)dr, (17) 
Jo 

where 

e^'\r) = l(-i)^'(i), 

for N^, i = 1, . . . ,d, some independent standard Poisson processes. Furthermore, we have 
recalled in Theorem 12.61 that the solution y to ([1]) is a continuous function of (a, i?,B^), 
considered respectively as elements of M*^, C7(M'^) and €2'^ {R'^^'^) for 1/3 < 7 < H. Thus 
our approximation theorem 11.11 can be easily deduced from the following result: 

Theorem 3.1. For any e > 0, let X^-^ = (Xff (z, j)) s,t>0;i,j=i,...,d be the natural Levy's 
area associated to X'^ , defined by 

X^^(z,j)= [\x^^-XndXi;^, (18) 

J s 

where the integral is understood in the usual Lebesgue-Stieltjes sense. Then, as e 0, 

where denotes the Levy area defined in Proposition \2.S\ and where the convergence in 
law holds m spaces (M"^) x Cl^ {M!^'"^) , for any ^ < H . 

The remainder of our work is devoted to the proof of Theorem 13.11 As usual in the 
context of weak convergence of stochastic processes, we divide the proof into the weak 
convergence for finite-dimensional distributions (Section [5]) and a tightness type result 
(Section m). 

Remark 3.2. A natural idea for the proof of Theorem 13.11 could be to use the methodology 
initiated by Kurtz and Protter in [18]. But the problem, here, is that the quantities we 
are dealing with are not "close enough" to a martingale. 

4. Tightness in Theorem 13.11 

From now, we write Cf (resp. C|^) instead of Cf(M'^) (resp. Cl^{W^'"^)). We first need 
a general tightness criterion in the Holder spaces and C^/. 

Lemma 4.1. Let S"^ denote the set of (x, x^) G C7 x C^^ verifying xq = and 

\/s,t> 0, V^, J = 1, . . . : x^,(2, j) = 4,{t,j) - - xlixl - xi). (20) 

Let fi such that < /i < 7. Then, any bounded subset ^ of S'^ is precompact in x Cg'^. 



Proof. Let (x'^jX^'") be a sequence of J2. By assumption, (x^jXq.'"") is bounded and 
equicontinuous. Then, Ascoli's theorem apphes and, at least along a subsequence, which 
may also be called (x",Xq.'"), it converges uniformly to (x,Xq.). Using ( l20l) . we obtain in 
fact that (x"',x^'") converges uniformly to (x,x^). Moreover, since we obviously have 

\\x\\n < liminf ||a;"'|L and ||x^||2u < liminf ||x^'"||2u, 

n— >oo n— >oo 

we deduce that (x,x^) G x Cg'^. Finally, we have 

— >0 and Ijx^ - x^'^Ha^, — ^ 0, 



owing to the fact that 

1 1 ^ 1 1 jj, 1 1 \ \ "y II 1 1 ^ 1 1 3j 1 1 '•y I j 1 2/ 1 1 'y ^ II ^ ^ 



and similarly: 

|2m < (||X^||27 + I|X^'"||27)||X- -X ' ||oo 



2, nil ^ /'ll„2|| I ||„2,n|| ^||„2 ,,2,n|| l-;^ 

□ 



We will use the last result in order to get a reasonable tightness criterion for our 
approximation processes and X^'^, by means of a slight elaboration of flUi Corollary 
6.1]: 

Proposition 4.2. Let and X^'^ he defined respectively by [n\ ) and fTE\) . If, for every 
f] > 0, there exists 'y > fJ^ and A < oo such that 

sup P[||X^||^ > A]<r] and sup P[||X2'^||2^ > A]<r], (21) 

0<£<1 0<e<l 

then (X%X2'=) is tight m x Cl^ . 

Proof. Recall the Prokhorov theorem relating precompactness of measures on a space 
to compactness of sets in the space. This result states that a family M of probability 
measures on the Borel sets of a complete separable metric space S is weakly precompact 
if and only if for every r} > Q there exists a compact set <Z S such that 

sup /i (S* \ Kr^) < 7]. 

Furthermore, it is readily checked that the couple (X^, X^'*^) satisfies the assumption 
( |20l) . which allows to apply Lemma [4.1[ Hence, combining this lemma with Prokhorov's 
theorem, our proposition is easily proved. 

□ 

Let us turn now to the main result of this subsection: 

Proposition 4.3. The sequence (X^, X^''^)e>o defined in Theorem \3.1\ is tight in xC2'^. 

Proof. Thanks to Proposition 14. 2[ we just have to prove that (X^, X^'^) verifies (pTl) . For 
an arbitrary G (0, 1), we will first deal with the relation 

sup POI^'IIt >A]<7], (22) 

0<e<l 



for A = large enough, and 1/3 < 7 < iJ. To this purpose, let us recall some basic 
facts about Sobolev spaces, for which we refer to P: for a G (0, 1) and p > I, the Sobolev 
space W"'^([0,T]; M") is induced by the semi-norm 



p 

a,p 



T rT 



Jo 



dsdt. 



(23) 



Then the Sobolev imbedding theorem states that, if ap > 1, then VV"'*'([0, T]; M.'^) is 
continuously imbedded in CJ{M.'^) for any 7 < a — 1/p, where the spaces CJ have been 
defined by relation ([8]), and in this case, we furthermore have that 



ll/|l7<C 



(24) 



for a positive constant c = Ca,p- Notice that, in both ([8]) and ( I23I) . the sup part of the 
usual Holder or Sobolev norm has been omitted, but can be recovered since we are dealing 
with fixed initial conditions. In order to prove (l22l) . it is thus sufficient to check that, for 
any p > 1 sufficiently large and a < H, the following bound holds true: 



sup E 

0<e<l 



T rT 



\X'{t)-X'{s) 

\f. — 



dsdt 



< Ma,p < 00. 



Invoking Lemma F6.lt we get, for any e > 0, any t > s > and any integer m > 1: 



< 



E [\X''\t) - X''\s)\'^"'] 



(25) 



(26) 



r)^^^ - (s + e 



rf--2)e''\r)dr 



2m 



{t + e-r)''-^e''\r)dr 



2m 



)m— 1 



< 



(2m)! 



m! 
2m-i 



(2m)! 



m! 



)m— 1 



< 



(2m)! 



m! 
2"!— 1 



(t 



\2H~l 



(s + e — r)^ ^^Ydr 



dr 



(2m)! 



m! 



) dr 



\t-S 



2mH 



{t - ry"-^dr 
by Lemma [2.71 



(27) 



(2^ 



Note that here, and in the remainder of the proof, C|.} denotes a generic constant de- 
pending only on the object(s) inside its argument, and which may take different values 
one formula to another one. From fl2Hl) . we deduce that fl2^ holds for any a < H and 
p large enough, from which ( l22l) is easily seen. Moreover, thanks to the classical Garsia- 
Rodemich-Rumsey lemma, see [12], for any e, 5, T > and i G {1, . . . , c?}, there exists a 
random variable G^'"''^'* such that, for any s,t G [0,T]: 



\X'^\t) 



X'''( 



< G 



T,5. 



,e,t I 



lH-5 



(29) 



Since the bound in fl28l) is independent of e, it is easily checked that, for any integer 
m > 1, any i E {I, . . . ,d} and any 6,T > {6 small enough), we have 



C2m,5 ■■= sup E{\G^^'^''' 
0<e<l 



Let us turn now to the tightness of (X.^''^)s>o- Recall first that Xg/(z, i) = 



Therefore, we deduce from 



that 



Cim,H 
02m 



t - S 



\imH 



(30) 



Assume now that i ^ ]■ We have, by applying successively (jSSl), Lemma EH] and (12^ : 



i?[|X?f(z,j 



|2ml 



< Cra.E 



+CmE 
+Cm,H E 



J s 

Jo 



,2H~1 



du 



u 



— {s + e — u)" ^^Ydu 



- Xl''\ {u + e- v)"-^du ) dv 



Vv 



This last expression can be trivially bounded by considering the case e = 0, and some 
elementary calculations then lead to the relation 



E[\X.^,f{i,j)r]<c^,H\t-s 



4mH—2mS 



(31) 



In order to conclude that verifies the second inequality in ( 1211) , let us recall the following 
inequality from [13]: let g G C2{V) for a given Banach space V; then, for any k > and 
p > 1 we have 



blU < c {U^+2/pA9) + W^gW-r) with U^-p{g) = (^j^ 





9st\ 


p 


\t 


— s 


7P 



dsdt 



i/p 



(32) 



By plugging inequality flHUjl - flHT]) . for 5 > small enough, into flH21) and by recalling that 
(5X^'^ = 6X'^ ® and inequality fl2^ . we obtain easily the second part of 0211) . 

□ 



5. FDD CONVERGENCE IN THEOREM 13.11 

This section is devoted to the second part of the proof of Theorem I3.H namely the 
convergence of finite dimensional distributions. Precisely, we shall prove the following: 

Proposition 5.1. Let (X^,X^'^) be the approximation process defined by and flE) . 
Then 

f.d.d. - lim(X^ X^'") = {B, B^), (33) 

e^O 

where f.d.d. — lim stands for the convergence in law of the finite dimensional distributions. 
Otherwise stated, for any k > 1 and any family {sj, ti] i < k,0 < Si < ti < T}, we have 

C - lim(Xf^ , X?;^^ ,...,Xl, X^f,^ ) = (B,, , Bl,^ , . . . , B?^,^ ) . (34) 



Proof. The proof is divided into several steps. 

(i) Reduction of the problem. For simplicity, we assume that the dimension d of -B is 2 
(the general case can be treated along the same lines, up to some cumbersome notations). 
For i = 1,2, e > and < n < t < T, let us consider 



J u 



and 

Y\u,t)= f{Bl-Bi:){v-uf~Uv. 



In this step, we shall prove that the fdd convergence fl33l) is a consequence of the 
following one: 

e''\u)du, I e'^\u)du, [ xl^'e''\u)du, 

Jo Jo 

Y'^^'{u,-)e''^{u)du, [ xl''e''\u)du, [ Y^''{u,-)e''\u)du 

Jo Jo 

(w\W\ [ BldWl, [ Y'iu,-)dWl [ BldWl [ Y\u,-)dW^) . (35) 
V Jo Jo Jo Jo / 

Indeed, assume for an instant that fl35p takes place. Then, approximating the kernel 
(t — ■)^~i/2 in by a sequence of step functions (along the same lines as in [TJ Proof of 
Theorem 1, p. 404]), it is easily checked that we also have: 



X''',X^'', / {■ + e - uf--^ Xl-' e''\u)du, 
Jo 

[ Y^''{u,-)e''\u)du, [ {■ + e - u)"-^ Xl'' e''\u)du, [ Y^^'{u, ■)9''^{u)du 
Jo Jo Jo 

{B\B\j'^{.-ur-'^BldWl 

[ Y'iu,-)dWl [ {--uf-'^BldWl [ Y\u,-)dW^). 
Jo Jo Jo / 

(36) 

In other words, we can add the deterministic kernel {■ + e — u)^~^ in the first, second, 
third and fifth components of ( l35i) without difficulty. Let us invoke now the forthcoming 
identity ( l53i) in Lemma [6.31 for s = 0, which allows easily to go from ( l36i) to: 

(Xl'^X2'^Xo'.•^(l,2),Xo'.'^(2,l)) (^B\B\J^B^dB\J^B'dB^y (37) 

Finally, in order to prove our claim flM|) from fl37|) . it is enough to observe that i) = 

{Xl'y/2 and 

X.?(z, j) = XoY(^, j) - Xo^f (z, j) - X^^{xr - Xr). 

(a) Simplification of the statement [3^) . For the simplicity of the exposition, we only prove 
( 135|) for a fixed t, instead of a vector (ti, . . . , tm)- It will be clear from our proof that the 



general case can be elaborated easily from this particular situation, up to some additional 
unpleasant notations. Precisely, we shall prove that, for any u := {ui, . . . ,Uq) G M^, we 
have lime^o'^e = E[exp{i{u, U))], where 6s := E[exp{i{u, U'^))], U'^ is defined by 

u' = ui[ e''\v)dv + u2 [ e'''^{v)dv + us [ xl''e''\v)dv 

Jo Jo Jo 

+ Uif Y^'%v,t)e''\v)dv + U5 f Xl''e''\v)dv + ue [ Y^''{v,t)e''\v)dv, 
Jo Jo Jo 

and 

U = uiWl + U2W^ + m [ BldW] 

Jo 

+ uJ Y\v,t)dW^ + u, [ BldW^ + u,[ Y\v,t)dW^. 
Jo Jo Jo 

In order to analyze the asymptotic behavior of (5^, let us first express f/^ as an integral 
with respect to 6*^'^ only. Indeed, Fubini's theorem easily yields 

xl''e'-\v)dv = [ due''\u) [ dve''\v){v + 6-u)^~-^, 



and the same kind of argument also gives 



t 

Y^'%v,t)9'^\v)dv 

t PW 



duO'^^iu) dw dve'''^{v){w-v)"~^{w + e-u)"-^ -(v + e-u)^-^) 



' u J u 

ft pu 

+ / due''^{u) I dw dve'''^{v){w-v)"-l{w + e-u)"-^. 

Jo Ju Jo 

Therefore, integrating first with respect to the randomness contained in 6^'^, one is allowed 
to write 4 = e^"^ /o ""'M'^'^) where, for / e L\[0,t]), we set 

and where the process is defined by: 

Zl ■= ui + u^Xl'^ + U4Y^''{u,t) + U5 [ {v + e-u)''-h''\v)dv 

J U 

rt rw 

+Ug dw dve''''^{v){w - v)^'i{{w + e - u)^~^ - {v + e - u)^"^) 

J u J u 

+M6 [ dw [ dve''\v){w-v)^-^w + e-u)^-^. (38) 



u)du 



Ju Jo 

Hence setting now, for / G L^([0,t]), 

f2, 



$(/) := E (^e'Iofi-)'iWi^ = exp {^-\ jj' 



we have obtained the decomposition 

Se = E (^$(Z)e'"^'^') + < + v'', 
where the process Z is given by 

Zu = ui + U3BI + UiY'^iu, t)+U5 I {v- u)"--^dW^ 

J U 

rt pw 

+U6 dw dW^ (w - v)^~^{{w - u)^~^ - (v - u)^) 

J U J U 

pt pu 

+U6 dw dW^ {w - v)"-^w - u)"-^ , ue[0,tl 

Ju Jo 

and with two remainders v^,v'^ defined as: 



The convergence of above is easily estabhshed: using again the same strategy than in 
[71 Proof of Theorem 1] (namely reducing the problem to a convergence of Kac-Stroock's 
process to white noise itself via an approximation of Liouville's kernel by step functions), 
one has that 

3^.2/ Law , 2^ 



z^ / e^'\u)du ^ iz,w,') 



e~>0 



Note that the convergence in law in the last equation holds in the space ^ x R, where 
^ = ^([0,t]) denotes the space of continuous function endowed with the uniform norm 
II ■ ||oo- In particular, it is readily checked that limg^o'^^e = 0. 

Now, it remains to prove that lim^^o = 0- To this aim, we notice that we can bound 
trivially |e*"^^* | by 1, and then apply the forthcoming Lemma [6^2] in order to deduce that 



\v^\ < E 



^2 4 U,|\ „2||2^||2, 



for any a G (0,1). Furthermore, it is well known that characteristic functions on a 
neighborhood of are sufficient to identify probability laws. Consequently, using Holder's 
inequality, we see that in order to get lime^o'i'e = 0, we are left to check that, for a given 
Mo > 0, 

sup E[\\Z,\\l] <oo, (39) 

0<e<l 

limE[0|^(e)] =0, limE[^|^(e)] =0, limE[^|^(e)] = 0, (40) 
sup E[e"'ll^^lli2] < M for all u < uq. (41) 

0<e<l 

We are now going to see that relations ( l39i) . ( 140|) and ( HTj) are satisfied. 



(Hi) Simplification of inequality Recall that Z"^ has been defined by f lHHj) . and 

decompose it as = ui + u^Uliu) + M4?7|(m) + u^U^iu) + uqUI{u) + uqU§{u), with 

U!{u) = Xl'^, UI{u) = Y''%u,t), Ul{u)= [\r + e-uf~'^e^'\r)dr 

J U 

pt PW 

Ul{u) = U6 dw dre'''^{r){w - r)^~^{w + e - u)"~^ - {r + e - u)^~^) 

J u J u 

pt pu 

Ul{u)= j dw j dr9'^^{r){w-r)^-^w + e-u)^~^ 

Ju Jo 

In order to obtain fl4ip . it is sufficient to check that there exists M > such that, for 
K, > small enough and z = 1, . . . , 5, we have 

sup E (e''-^oUH^fdu\ ^42) 

0<£<1 ^ ^ 

Moreover, observe that t/f can always be written under the form 

Ul{u)= I Vi{u,r,e)e''\r)dr, (43) 
Jo 

for a deterministic function Vi{u,r,e), and it is thus enough to check that 

Ci := sup sup / V-^ {u, r, €)dr < oo. (44) 
«e[o,T] o<£<iJo 

Indeed, using Lemma 16.11 below, we can write 

m=0 L ^"^0 ^ - 



-m=0 m=0 

where we have used the bound (m/S)™" < m! < in the last inequality, so that the 
desired conclusion follows for n > small enough. 

(iv) Proof of ^4\ l- We shall treat separately the cases for i = 1, . . . , 5. During all the 
computations below, C > will denote a constant depending only on H and T, which 
can differ from one line to another. 

a) Case i = 1. We have X^'^ = Vi{u,r,e)9'^''^{r)dr with 



\/i(u,r,e) = l[o,„](r)(n + £-r)^ 2 



T „,2H 

Hi 



Since 

V^{u,r,e)dr= I {u + e - rf^'^dr < I {u - rf"-^dr = ^ < C, 
Jo Jo 2ii 

we have that ( l44l) takes place for i = 1. 



b) Case i = 2. We have Y'^''^{u,t) = V2{u,r,e)6'^''^{r)dr, with 



V2{u,r,e) = l[o,n](r) / [{w + e - r)^ - {u + e - r)^ ^^{w-u)^ ^dyj 

J u 

+ l[u,t]{r) / (w + e - r)""^w - u)^~'^dw. 

J r 

Then V2^(u, r, e)dr = A2a{u, e) + A2^2{u, e), where 

A2^i{u,e) = J [{w + e — r)^~^ — {u + e — r)^~^){w — u)^~^dw^ dr, 

A2,2{u,€) = J i^j {w + e - r)^~^{w - u)^~^diu^ dr. 

For any j3 G (0, 1) and > n > r > 0, we can write, for some w* G {u + e,w + e): 
\{w + e - r)^"^ - {u + 6 - r)^~^ | 

C\w - u w - 



— r|^2 ^^1^ \\w + £ — r\2 ^ \u + e — r\2 ^ I |m — r|2+^ ^ 

Then, choosing P = ^ — H + 6 (with 6 > small enough), we can write 



t . N 2 



AiK^) < CI ^-Iz^xf/ ^^^1 <C, 



where we have used the fact that 2 — 4H < 1 whenever H > 1/4. On the other hand, 
using that, for a G (0, 1), > 1 and A > 0, we have 

dy I r dy C 



we obtain that 

A2,2{u,e) < C fif- dr 

c) Case 2 = 3. We have 

{r + e-u)"~h'''^{r)dr= / V3{u,r,e)e''^{r)dr, 

Jo 

with V3{u, r, e) = l[u,t]{r){r+e — u)^"^, so that the desired conclusion follows immediately 
since 



y3'(M,r,e)rfr = / {r + e - uy'-'dr < / (r - uy'-^dr = ^i/ - 



d) Case i = 4. We can write 



u J u 



dw I dr{w - r)^~i [{w + e - u)"~^ - (r + e - u)"~^)9''''^{r) 



as V4(m, r, £:)6'^'^(r)(ir, with 



/"* 3 1 1 

V4{u,r,e) = l[u,t]{r) / {w - r)^^^ [{w + e - u)"^^ - (^r + e - u)^~2^dw. 

J r 

Then, according to the computations already performed for the analysis of A2^i above, we 
obtain, for 6 > small enough. 



T rt 

Vi{u,r,e)dr < C 



1 



_ ,,|2-4H+25 



r — M 



dw 



\w — r 



i-s 



dr < C. 



e) Case i = 5. We have 

dw / dr{w — r)^~'2(^itj + e — u)^~'2 0'^^'^(^r) 
Ju Jo 

with 

H- " 



V^{u,r,e)d^''^{r)dr. 



l^5(u,r,£) = l[o,„](r) (w-r) ^(^w + e-u) ^dw. 

J u 

Since \w — r\^~^ > \w — u\^~^~^^\u — r\^~^ for r < < w,we get (for 5 > small enough) 
that 



w — r) '^[w + e — u) 2 



< 



C 



\w — r\2 \w — u\'i 



< 



c 



m — r 2 "\w — u\'2 



Hence, invoking again the fact that if > 1/4, we end up with 

dr I /"* dw 



1*1 i*u 

/ V^{u,r,e)dr <C 
Jo Jo 



u — r 



1-25 



W — M 2 



< C. 



(v) Proof of In the previous step, we have shown in particular that, for any i 

1, ...,5, we have supQ<g<;^ E[\Ui{u)\P~\du < oo for all p > 1, which implies 



sup / ^'[iZ.fJPjc/M < oo, for all p> 1. 

0<e<l Jo 



On the other hand, a simple application of Schwarz inequality yields 



iZD'e-'^du 







< 



Jo 



and the same kind of argument also gives 

E[^Ue)\=E 



dx r dy{Zl)\Zlfe 
Jo 



<'-E 

- 2 



dx / dy{Z'Xe 



^0 



2' 














2 





<Ce' I E[{Z-f]du. 
Jo 

Finally, we have 



E[^Ue)\=E 



e\\Z^ 



\ 1/2' 

{Zlfdu 

/ 



1/2 



<2e'Em\\l.)^2 E[{Zlf]du<2e'E{\\Z^\\l^ + 2e'l^ i^j^ E[{ZlY]du 
and the proof of (HUjl follows immediately by putting all these facts together. 



(vi) Proof of For all a< (3-^, the Sobolev inequality (El]) yields \\Z^\\a < C\\Z^\\p^p, 
where ||/||/3,p has been defined by ( 1231) . Moreover, recall from (l38l) that Z^ has the form 

G(s,t,r)^"'2(r)c/r 



for some (j(s,t, ■) G L^([0,T]). Hence, using the definition of 6'^'^, we can write, for any 
even integer p > 2, 

E\Zt-Zl\^ 

= e-^ f G(s,t,ri)---G(s,t,rp)E[(-l)^('^)+-+^(^)]dri---rfrp 

J[0,T]P 

_ f _2{ri-r2) _2(rp_i-rp) 

= pie ^ G{s,t,ri) ■ ■ ■G{s,t,rp)e ^ ■■■e 7^ l{r^>...>r^}dri ■ ■ ■ drp 

J[0,T]P 

£ " I G{s,t,ri)G{s,t,r2}e 7^1^^._^^^^ydridr2 

'[0,T]2 



(p/2)! 



(p/2). V- G(M.rOG(..Vr.)S[(-l)-'-'+«'-']*,*. 



2\P/2 



2Hp/2)\ 

In particular, we see that, in order to achieve the proof of ( l39ll . it is enough to check that 

E\ZI - Z'^\^ < C\t - sf~^ (45) 

for some 6 > (small enough). Actually, we shall use again the decomposition of in 
terms of the f/j's, which means that it is sufficient to prove E'If/f (m) — f/f (f )p < C\u—v\^^^ 



for i = 1, 5. But it is easily seen that 



E 



T 

{Viiu,r,e)-Viiv,r,e))e'^\r)dr 



< 



T 



{yi{u,r,e) - Vi{v,r,e)Ydr, 



where Vt is defined by ( 143|) . and are specified at step (iii). It is thus enough for our 
purposes to show that J^{Vi{u,r,e) — Vi{v,r,e))'^dr < C\u — v\^~^ for < m < u < t, 
where t G [0, T], and we will now focus on this last inequality, for i = 1, 5. 

a) Case i = 1. We have 

(yi{u,r,e) — Vi{v,r,e)) dr 



[{u + e - r)^ 2 _ [I, 



e — r 



+ e - r)^~^Ydr + j {u + 

J V 

pv pu 

< / ((m - r)^^^ - (t; - r)^~^)^rfr + / {u-rY"~^dr 

Jo Jv 

<c{uf^ by Lemma O 
b) Case i = 2. We can decompose (V2(u,r, e) — V2{v,r,e)Ydr into 

[ {V2iu, r, e) - V2{v, r, e)fdr = B2,i + ^2,2 + ^2,3 + ^2,4 
Jo 



,2H-1 



dr 



with 



-82,1 



\Ju 



u / pt 



[{w + e — r)^ 2— (w + e — r)^ ^2){w — u)^ ^dw 

t ^ ^ ^ x2 

[{w + e — r)^~2 _ -I- 5 — r)^~2^(^itj _ v)^~^dw ) dr 



B 



2,2 



V \J u 



{^{w + e — r)^ 2 — [u + e — r)^ 2^{w — u)^ 2dw\ dr 



and 



B. 



2,3 



{w + e — r)^ ^{w — v)^ ^dwj dr 
{w + e — r)^~^ [{w — u)^~^ — {w — v)^~i^dw j dr. 



i?2,4 = 

We will now study those terms separately, starting from -82,1: this first term can be 
further decomposed as i?2,i < c(i?2,i,i + -82,1,2 + -82,1,3), with 

-82,1,1 

[{w + e — r)^ 2 — [u + e — r)^ 2) x {{w — u)^ ^ — (w — v)^ ■2)dw j dr 
\ Ju j 



and 



( / [{w + e — r)^~2 — (^v -\- e — r)^~2^{w — v)^~2dw\ dr. 

Moreover, for 5 > small enough, and by choosing (3i = H — 5 and P2 = y ~ 
write 



B2X1 < C I \ I : ^7^\^-uf'- 



\Ju It; -r|^i(i-^) 



2 

\u — vl^'^ -r, -dw 1 dr 



since | — — S < 1 whenever H > 1/3. The term -82,1,2 can be bounded in the same 
way: for 6 > small enough, and by choosing again jS = H — 6, we can write 

^2,1,2 < C [ I [ La m l^-^l^i ,n mn m , "^l?^^) 

< C\u - vf-'' r \— ( f ^-^^TT-.dw] dr < C\u - vf-''. 

Finally, for 6 > small enough, and still with f3 = H — 6, we can write 

^ '■^ 1 , ,^ 1 1 ' ^ 



^2 13 < C \ — -\w-vf -J 3 — dw] dr 



fv 1 ■ -2 

< C I :-— I / -dw 1 dr 



Q \V — r\ \J^, y; _ ^ 2 



< C|«-.r"--y^ j-—^dr<C\u-v 

By gathering the bounds we have obtained for i?2,i,i, -^2,1,2, -^2,1,3, we now easily get i?2,i < 
C|m — v\^^^ for any arbitrarily small constant 6 > 0. 



Let us turn now to the bound concerning -82,2- for j3 = ^ — H + 6 with 6 > small 
enough, we can write, by using similar arguments to those used in the study of -82,1,3, 

^2,2 < C I [ I } ^-^dw\ dr 

< c T: ^^^( f -. ^-^dw\ dr<C\u-v\"-^'. 



In a similar manner we obtain, still for 5 > small enough, 
i?2,3 < C r ( f \^dw\ dr 

t 1 



\H-25 



^ ^ I \^_r\2-AH+28 [J^ ) dr<C\u-v 

Finally, by choosing /3 = -j — 25 with S > small enough, the same kind of arguments 
allow to write 

^2.4 < C I ( [ L^\u-v\^- -^^^dw] dr 



In conclusion, putting together the bounds we have obtained for all the terms 5, we 
end up with the announced claim (V2(M,r, e) — V2(f,r, e)) dr < C\u — v\^~^ for any 
< S < H, as soon a.s H > 1/3. 

c) Case i = 3. This case can be treated along the same lines as the previous ones: we 
have 

{V3{u,r,e) -V3{v,r,e)) dr 

(l[„,t](r)(r + e- n)^"^ - l[^,t](r)(r + e- v)"-^Ydr 
<C( [ {r-v)^^-^dr+ [ {{r-uf-^-{r-vf-^YdrY 



Now it can be shown, by using the same reasoning as in the proof of Lemma 12.71 that 
Jq [Vsiu, r, e) — Vslv, r, e)Ydr < C\u — v\^^ . 

d) Case i = 4. Let us decompose again the integral under consideration into: 

/ {V4{u,r,e) - V4{v,r,e)Ydr = B^^i + ^4,2, 
^0 



with 



^4,2 = j (^j^{w -r)^~i{{w + e - u)"~^ - {r + e - u)"~^ 

1 1 \ ^ 

— {w + e — v)^~'2 + {r + e — v)^~'i)dw\ dr. 

The first of these terms can now be treated by using similar arguments to those used in 
the study of 53^2: ioi (5 = \ — H + 5 (with 5 > small enough), we can write 

Bii < C I I / „A w-rf } ^ dw\ dr 



V \ J r 



t 1 



< C / TTTTTT^ / 1 TT-.dw dr<C\u-v 



\H-2& 



In order to estimate -84,2, we introduce an interpolated point (x*, y*) G {r + e,w + e) x 
(t;, u) which allows to write, for u < r, 

\(w + e — u)^~^ — ir + e — u)^~^ — (w + e — v)^~^ + (r + e — v)^~^ I 



< [-- H ] {-- H ] \x* -y*\"~^\w -r\\u-v\ < C 



\w — r\\u — v\ 



I I £_ W 

r — n 2 



Therefore, with /5 = -|- — 5 and 5 > such that | — |if + 5 < 1 (recall again that this is 
possible as soon as if > i), we get 

2 



t / rt 



u \ J r 



1 1 ,, 1 



B42 < C \ 3 -5 \\w-rf\u-v\'^ -J dw\ dr 



t ■ ■ ^2 



I I 3 3 ZJiS 

lyj — 2 2 



2 



e) Case i = 5. Our last decomposition can be written (V5{u,r,e) — V5{v,r,e)Ydr < 
C(fi5,i + 55,2 + 55,3), with 

i?5,i = J (tf — r)'^~t(-u; + £ — M)^~^du;^ dr 

-65,2 = / ( / (tf - r)'^"^ ((u; + e - n)^~^ _ (y; + e _ t>)^~5)c?y; j c?r 



^0 V-'u 

i?5,3 = y (u; — r)^"^(w; + e — 1 (ir. 



■u / f u \ 2 

3 , ^ cj 1 



Then, choosing 6 > small enough, and thanks to the fact that 
for all < r < n < < t, we have 

<C r 1— — ( f ^—r-,dw\ ' dr < C\U - < ^1 _ |H-25^ 

On the other hand, for /5 = y — 25 with 6 > small enough, we can write, using similar 
arguments to those used in the study of -82,3? 

\v — r\ \^Ju \w — u\2 2-" 
Finally, still for 6 > small enough, invoking again inequality fHUl) . we have 

B,,3 < C r( T- \^dw] dr 

Jo \Jv \w — r\'2 " \w — v\2 J 

2 

i / / i , \ . _ „, 'H-2& 



< C , ,12-3^-2. / ] ^T^^^ dr<C\u-v 

Jo 1^' — '^l \Jv \W—V\ 2+" ' 

Our estimates on -85,1, -85,2, -85,3 yield trivially {y^{u, r, e) — V5{v, r, e)Ydr < C\u—v\^~^ 
for S arbitrarily close to 0, and this last bound finishes the proof of our Proposition 15.11 

□ 



6. Some technical lemmas 



This section collect the technical results that have been used throughout the proof of 
Theorem 13.11 The first lemma aims at giving some estimates concerning the Kac-Stroock 
kernel which can be seen as a elaboration of the ones contained in Delgado and Jolis 
[3, Lemma 2]. Notice however that these latter results are not sharp enough for our 
purposes, which forced us to a refinement. 



Lemma 6.1. Let m G N, /, /i, 

2m ftrp 



,f2m e L2([0,T]), k e {1,2} ande> 0. We have: 



E 



n 



f,j{r)6^'^{r)dr 







<^ll/l|U2...||/: 



and 



E 



f{r)9''^{r)dr 



2m+r 



2'"m! 



(2m + 1) 
2™+im! 



2m||L2, 



2m 
L2 ) 



(47) 



(4^ 



where 



Proof. For m e N, e > and fi, . . . , f2m & L'^{[0,T]), let us denote 



2m) 



E 



2m. „y 

n / fAr)0''\r)dr 

.7=1 



We will need to introduce some operations on the set of permutations (in the sequel, 
&k stands for the set of permutations on {1, . . . , A;}): when r G &2m and cr e &m, we 
note cr -kT the element of ©2m defined by 

(a*r)(2j-l)=r(2a(j)-l) and (a ^ r)(2j) = r(2a(j)). 

Remark that we have id t*: r = r and cr' ^ (cr ^ r) = {(Jcr') * r, so * : &rn x '52m &2m 
defines a (right) group action of &m on ©2m- For any r G 62m, the orbit of r has exactly 
m! elements. Consequently, the set of the orbits under the group action -k has l^ZIlll 
elements and we have, by denoting Tj one particular element of the orbit Oj = o(ri) G ^: 
for ri, . . . ,r2m G [0, 1], 

l{Vi^J, nT^r-j} = l{»-r(l)>--->?-r(2m)} ~ {''Vdl > ■ ■ ■ >?-r (2m) } 

TG©2m OiG^rSOi 

(2™)!/m! m 

^ E ni{^-.0)-i>^-.0)}- (49) 

For the reader who might not be completely convinced by this inequality, let us illustrate 
it by an example: when m = 2 and Tj = id G ©4, we have Oj = o{Ti) = {id, (13)(24)} and 
we have used 

^ ] -'-{^■r(l)>»"T(2)>»"T(3)>'"T(4)} ~ I{ri>r2>r3>r4} + I{r3>r4>ri>r2} ^ I{ri>r2}l{r3>r4} 

2 

= n l{'-2r,U)-l>'-2riU)}- 

Let us apply now inequality ( l49l) . We introduce first a notation which will prevail until 

the end of the article: for e > and r G M+, we set Qeir) := je^ . Notice then that, 
for any e > 0: 

I^2m(/l5 • • • 5 /2m) I 
1 



< 



|/l(ri)|...|/2m(r2m)| E (-l)MA^(^) 



2m 



X^X^TS^ / lK(i)>...>r,(2m)}l/l('^l)l ■ ■ ■ l/2m('^2m)| 

/ m 

^{rr{2i-i) - rr{2i)) ) ciri ■ ■ ■ cir 



1=1 



and thus, according to fH9|) . we obtain 



^2m(/l5 ■ ■ ■ 5 hm 

(2»n)!/m! ^ 

< 

(2m)! 



< 



'/••-■ III, n 

5^ n / l{ri>r2}|/2r,(i)-i(ri)||/2r.(j)(r2)|Qe(ri -r2)c/rirfr2 

II/1IU2 . . . ||/2m||L2, 



2"*m! 

the last inequahty coming from 

I{ri>r2}l/fc(n)||/^(r2)| Qe{ri - r2)drxdr2 

[0,T]2 



|/£(r2)r / Q.(ri-r2)t/riUr2 



(50) 



This finishes the proof of (H7I) . so let us now concentrate on (HHl) . For m G N, we have 



E 



rp ^ 2m+l' 



f{r)e''\r)dr 



< 



2m+l 



n i/(^')i F 



1=1 



dri . . .dv' 



2m.+l 



2m + 1 

c-2m+l 



|/(s)|rfs 



T 



2m 

ni/w 



(iri . . .dv' 



2m 



< (2m+l)AUI/l,---,l/l) / 1/(^)1-6-1- ci^. 







1—2" 1 

Since for s > e, we have that — e ^ < we get that 



\fis)\-e-"^ds 



< { 1^ imi'ds 



and (jlH]) follows easily. 



\f^s)\-e-"^ds + e ! \f{s)\\e-"^ds 



(51) 



Qe{2s)ds 



1/(^)1^^ <o 1/(^)1'^^ +0II/IU2 



□ 



The following lemma aims at measuring the distance between the laws of the stochastic 
integrals /J" f{r)Q^'^{r')dr and f{r)dW^, whenever / is a given (deterministic) function: 



Lemma 6.2. Let f G C"([0,T]) for a given a G (0,1), k G {1,2} and e > 0. For any 
u eM., we have: 



< 



^ '^a \\J \\a 



U\ 



L2U^ + (pf{e)— + i^fie)— + ipf{e) ^ 



with Ca = /q°° "^^dx and 



e 2 



(52) 



dx / dynx)f\y)e-^ 

^0 

Proof. The proof is divided into two steps. 

1. First step: control of the imaginary part. We can write, thanks to (148!) : 



< 



m=0 



2m+l 



(2m + 1)! 



7^ ^ 2mH-l' 

£,fc 



f{r)e''''{r)dr 



\u\ 



" "-"Ira 



2. Second step: control of the real part. This step is more technical, and we will mainly 
get a bound on the quantity Lm,e defined by: 



f'^{Sm)dSr^ 



In order to express this quantity in a suitable way for estimations, notice that e "^^ds = 
|. We can thus insert this term artificially in the multiple integrals involved in the 
computations of £'[e*"-^o /{'')'^'^r ]. This gives: 



1 



(2m) 



AL(/, •••,/) 



f (ri)dri / e-^'~'dr2... I f\r2m-i)dr2m-i I e-'^^'-dr 



r2m-3 



-2r2„ 



2m 



By a telescoping sum argument, we can now write L^^^ as a sum of m terms, whose 
prototype is given by M„i^£ = M^ ,, + M^^^ - M^^^, with 



m,e 



ri 



\f{ri)\dri / \f{r2)\Qe{ri-r2)dr2... 



\f{r2ra-\)\dr2ra-\ / | /('^2m) " /('"2m-l) | (^^2 m— 1 ^2m 

where ^ is defined by 



f{ri)dri / f{r2)Qe{'ri - r2)dr2 . . . 



r2m.-2 



f {r2m.-i)dr2 



m — 1 



Qe{r2m) dr 



2m y 



and where 



Jo Jo Jr0.rr,.-7, Jo 



We will now bound those three terms separately: invoking first fl50|) . we get 

X \fir2m-l)\dr2m-l\\f\\a ^2™ " ^^2™-! r<5e('^2m-l - ^"2™) C?'"2m 

7o Jo 
On the other hand, ( l50i) and ( ISTl) also yield: 

2 ^ 



< ^^^^-^ (^y^ \fin)\dn I \f{r2)mn-r2)dr2 







X llf{r2m-l)Qs{r2m-l)dr2m-l < _\ ) ,3^ II / II V/ (g) 

Finally, ^ can be bounded in the following way: observe that 



1 



JO 



X \f{r2m-3)\dr2m~3 / | /(r2m-2) |Q£(?^2m~3 " '"2m-2) C??^2m-2 / / (?^2m-l )c?r2m,-l , 

^ Jo Jo Jr2m-2 



which can also be written as 



Mi,< 7 %; r||/|P,T"^M^'l, 



with 



= / dr2m.^lf{r2m~l) / C?'^2m-2 / C?^2m-3 1 /('^2m-2) 1 1 /('^2m-3) I 

Jo Jo JT2m-l 

2(r2 

771 — 3-'"2m-l) 2{r2 

e e 

X 



It is now readily checked that 
< [f{r2m-i)dr2m-i 



r2m-l 

dr 



2m-2 



2m.-3 



''2m -1 



X 



2(r2 

77T — 3 ''"2m — l) 2{r2 

2/'^_ „ 



< 7; fir, m—l 



„ 2('-2m-l-'-2m-2) 
2^», N„ 



dr2ni-2f {r2m-2)e 



m — 3 ^2m— 1) 

+ / dr2m-3f (r-2m-3)e 



Consequently, 



< i 

'"'^ - (m-2)!2'^ 



Our proof is now easily finished: plug our estimates on M^^, and into the 
definition of M^ ^, and then in the definition L^ g. This yields 



„,2m 

< E 



m- 



(2m)! 



Ak(/,---,/)- 



(2m)! 
2™m! 



f'^{r)dr 



< 



which is our claim. 



^^"cJ|/|U||/|U2n2 + </>Ke)^ + ^^(e)^ 



,.2„,||2^ 



e 2 



□ 



The following lemma gives an alternative form for X^'^ and B^: 
Lemma 6.3. Fix i,j G {1, . . . , d}, and t > s > 0. For all e > 0, we have 

x?r(^,j) 

= [ {Xl^' -Xi^'){t + e-u)^'-^9''\u)du- [ {Xl^' - Xi^'){s + e - uf-h''\u)d 
Jo Jo 

-an! dve''\v) [ du{Xl'' - Xi''){u + e - v)""-! , 



u 



(53) 



where we have set an = 1/2 — H . In the limit e ^ 0, we also have 



{Bi-Bl)it-u)''~'^dW: 



-an I dWl / du{Bi-Bl){ 

JvVs 



{Bi - Bi) [it - «)^"^ - (. - «)^-i] dW: 



u — V) 



(54) 



Proof. For any e > 0, the process X^'* is different iable, and according to ffTTl) . we have 

Jo 

Recall also that we have set SX^f = Xj'^ — X^''^ for any s, t G [0, T]. This allows to write: 
^'fihj)= f 5XitdXl:^ = e"-i f 5Xite^^\u)du 

J s J s 

-an du6Xi^^ dv{u + e - vf-k'^\v). (55) 



Moreover, an elementary application of Fubini's theorem yields: 

du 5Xi^^ / dv{u + e- v)^-^^e^^\v) = / dv9''\v) / du SX^^ (u + e- v)"-l 
Jo Jo J syv 

= [ dve''\v) 5X1:^ [ dr{r + e-uf-l+ [ dve''\v) [ dr 6Xif{r + e - vf-^ 

Jo J sWv Jo J syv 

Integrating the kernel (r + e — u)^~^ , and plugging the last identity into ( l55i) . we obtain 
the desired relation (l53l) . 

To get formula fl54|) for B^^(z,j), it suffices to observe that 

BUz,j) = -hmjjBi- BDdBr 

with i?^'^ = Jq{u + e — v)^~^dW^, and then to mimick the computations allowing us to 
write (13^ just above. Details are left to the reader (see also the proof of [31 Lemma 3]). 

□ 

Finally, the following lemma gives an estimate for the variance of B^j(i,j) which is 
useful in the proof of Proposition 12.81 

Lemma 6.4. There exists a constant c > 0, depending only on H , such that B^j(z, j) 
< c\t — s|^^ for all t > s > and i,j G {1, . . . ,d}. 

Proof. The case where i = j is immediate by Lemma 12.71 so we only concentrate on the 
case where i ^ j. Using formula ( |5^ . we see that we have to bound the three following 
terms: 

Ai := E\Bi- BlWt-uf"-^du 

A2 := E\Bl-Bl\^ (^{t-u)"~-^ -{s-uf^-^y du 



A, := f du{Bi-Bi){u-v)^' 

Jo JvVs 



2 

dv. 



Throughout the proof, c will denote a generic constant (depending only on H, T) whose 
value can change from one line to another. Owing to the fact that E\Bl — Bi\'^ < c|n— sp^. 



t pt 



2 

du 



see Lemma 12.7^ we can write 

< c [\u - sf'it - uf^'-^du < c{t - sf"" / it - uf"^^ = c{t-s 
We also get 

A2 < c j\s-uf'' ((t-u)''--^-{s-uf-^) 
= c j^u^^ [{t- s + u)"~~^ - u^^^y du 
= cit - sY"" ((1 + uf~-^ - n^-^) ' du 

+ du = cit-s) 

the last integral being finite since H < \. Finally, we have 



E 



f du{Bi-Bi)iu-vf- 
du f dwE[{Bi-Bi){Bi,-Bi)]{u 

JvVs 



vVs JvVs 

< cf du f dw{u-vY"-^{w -vf"-^ = c( I {u-vf^"^du 
SO that 

[{t - _ (s - 'dv + c J' (^j\u - vf^'-'idv^ ' dv 



^3 < c 

< cit - s)^^ j [(1 + - ^2^^-^] 'dv + c I {t- vY^'-^dv 

This finishes the proof of the lemma. 
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